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We present a theoretical framework allowing to properly address the nature of surface-like eigen- 
modes in a hypersonic surface phononic crystal, a composite structure made of periodic metal 
stripes of nanometer size and periodicity of 1 /im, deposited over a semi-infinite silicon substrate. In 
surface-based phononic crystals there is no distinction between the eigenmodes of the periodically 
nanostructured overlayer and the surface acoustic modes of the semi-infinite substrate, the solution 
of the elastic equation being a pseudo-surface acoustic wave partially localized on the nanostruc- 
tures and radiating energy into the bulk. This problem is particularly severe in the hypersonic 
frequency range, where semi-infinite substrate's surface acoustic modes strongly couple to the pe- 
riodic overlayer, thus preventing any perturbative approach. We solve the problem introducing a 
surface- likeness coefficient as a tool allowing to find pseudo-surface acoustic waves and to calcu- 
late their line shapes. Having accessed the pseudo-surface modes of the composite structure, the 
same theoretical frame allows reporting on the gap opening in the now well-defined pseudo-SAW 
frequency spectrum. We show how the filling fraction, mass loading and geometric factors affect 
both the frequency gap, and how the mechanical energy is scattered out of the surface waveguiding 
modes. 

PACS numbers: 62.25.-g, 68.35.Iv, 43.35.+d, 43.35.Pt 



I. INTRODUCTION 



The idea of creating solids with artificial modula- 
tion in their physical parameters^ has proved fruitful 
both in the frame of coherent electronic transport and 
photonics, with the advent of semiconductor superlat- 
tices and photonic crystals, respectively. The anal- 
ogy between photons and phonons suggested considering 
phononic crystals^ periodic elastic composites of two or 
more vibrating materials, the interest being triggered by 
the possibility of achieving a frequency gap in the elastic 
modes dispersion relations and the wealth of applications 
stemming from tailoring the gap itself. 

Extensive investigation of acoustic band structures of 
periodic elastic composite has been carried out on two- 
dimensional phononic crystals both infinite^ and surface- 
terminated ^^i^ The archetypal 2D structure consists of 
parallel rods embedded in an elastic background. Rods 
can be actual or virtual ones (drilled holes). In the 
real case of a surface-terminated phononic crystal, an 
interesting class of acoustic modes, addressed as sur- 
face acoustic waves (SAWs) - arising from breaking of 
the translational symmetry when passing from an infi- 
nite to a semi-infinite medium - propagate confined to 
the elastic medium surface, the penetration depth be- 
ing of the order of their spatial period. Recent works 
elucidated some of the characteristics of surface acoustic 
waves in 2D phononic crystals of different geometrical 
configurations J^J^ 

The present paper focuses on pseudo-surface acoustic 
modes, reminiscent of surface acoustic waves, in surface- 
based phononic crystals. These structures are obtained 



patterning the surface of a substrate material either by 
metal deposition or substrate etching. The influence 
of a periodically structured overlayer on the acoustic- 
field eigenmodes of a thick or thin homogeneous slab 
has been studied by many research groups and effects 
such as band folding, ^^^^^ mode leakage, opening of fre- 
quency gaps^ — and interaction between slab and over- 
layer laode^^— have been discussed. Contrary to 2D 
phononic crystals, for surface acoustic modes in surface- 
based phononic crystals there is no distinction between 
the eigenmodes of the nanostructures and SAWs, the 
solution of the elastic equation being a pseudo-surface 
acoustic wave partially localized on the nanostructures 
and radiating energy into the bulk. This issue is particu- 
larly severe in the hypersonic frequency range, where the 
semi-infinite substrate's surface acoustic modes strongly 
couple to the periodic overlayer, thus preventing any per- 
turbative approach. In this physical scenario, we present 
a theoretical framework allowing to properly address the 
nature of surface-like eigenmodes in a hypersonic sur- 
face phononic crystal. We solve the problem introducing 
a surface-likeness coefficient as a tool allowing to find 
pseudo-surface acoustic waves and to calculate their line 
shapes. Having accessed the pseudo-surface modes of 
the composite structure, the same theoretical frame al- 
lows reporting on the gap opening in the now well-defined 
pseudo-SAW frequency spectrum. We show how the fill- 
ing fraction, mass loading and geometric factors affect 
both the frequency gap, and how the mechanical energy 
is scattered out of the surface waveguiding modes. 

Our results are of interest also from an applica- 
tive stand-point. The quest for hypersonic acous- 
tic waveguides, sources of ultrafast coherent acoustic 



waves and nano-opto-acoustic transducers in general op- 
erating at ever higher frequencies, requires patterning pe- 
riodic structures of ever decreasing periodicities. In this 
frame, surface-based phononic crystals provide a techno- 
logical advantage over 2D phononic ones, the processing 
technique involving the surface only, being easily scalable 
below 100 nm and suitable for high frequency transducers 
technology. The theoretical tool here presented, together 
with the comprehension of how the construction param- 
eters affect the frequency gap and surface waveguiding 
mode, will be a valuable tool for inspecting the pseudo- 
acoustic modes relevant for applications in view of device 
engineering beyond a trial-and-error approach. 

The present work is organized as follows: in Sec. ITTlwe 
outline and solve the full mechanical problem, calculating 
the vibrational normal modes of the composite system. 
The presence of periodic nickel stripes modifies the prop- 
erties of surface waves. ^^^^ The SAW is no longer an 
eigenmode of the elastic eigenvalue problem. Proper as- 
sessment of a pseudo-SAW in a surface-based phononic 
crystal is achieved in Sec. Illli We introduce a SAW- 
likeness coefficient permitting to discriminate pseudo- 
SAWs among the entire set of eigenmodes and to cal- 
culate their line-profiles. Investigation of the gap open- 
ing in a surface-based phononic crystal as a function of 
the filling fraction p, beyond a perturbative approach, is 
presented in Sec. |TVl We show how the frequency gap 
depends on the set of geometric factors {p, h} and mass 
loading {m, p} - p and h being the nickel stripe's mass, 
density and height, respectively. In Sec.|Vlwe analyze the 
energy distribution in the system as a function of p and 
relate it to the scattering of SAWs into the bulk. 



II. MECHANICAL MODEL 

We assume the system to be an elastic continuum com- 
posed of a periodic metallic grating deposited on a silicon 
substrate. In Fig. [1] we consider the general configura- 
tion where elastic isotropic nickel stripes of width d and 
height h = 50 nm are deposited on a crystalline silicon 
substrate. The grating has period A = 1 /im, thus ensur- 
ing pseudo-SAWs in the hypersonic range. Calculations 
are performed increasing (i, in order to cover the entire 
filling fraction p range (0, 1) - where p = d/X - and ex- 
plore different regimes, from a perturbative one (p <C 1) 
to substrate full coverage {p = 1). The acoustic equa- 
tion of motion governing the displacement u(r, t) of the 
composite system is 

dj [Cijmn (r) dnUm] = p (r) Ui , (1) 

where p (r) and Cijmn (r) are the position dependent mass 
density and elastic stiffness tensor, respectively, and the 
summation convention over repeated indices is assumed. 
For an harmonic time dependence e*'^^ + c.c., Eq.[T] reads 
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FIG. 1: Schematic diagram of the surface phononic crystal. 
Isotropic nickel stripes of width d and height h are deposited 
on a crystalline silicon substrate. The grating has period A. 
Propagation direction is along the x axis, surface normal to 
the bulk is along the y axis. 



We solve the eigenvalue problem via the finite elements 
method. As shown in Fig. [2fa), the model consists of 
a two dimensional silicon rectangular cell with a nickel 
stripe on top. The material properties for the silicon sub- 
strate and the nickel stripes, that enter in the expression 
of the elastic stiffness tensor, are reported in Table HI The 
silicon substrate's crystalline orientation is accounted for 
in the expression of the elastic stiffness tensor and, in the 
present case, the x axis is taken along the Si(lOO) crys- 
talline direction. To reproduce the entire nanostructured 
composite from the single unit cell, the displacements ui 
and U2, calculated respectively on side 1 and 2 of the 
cell (see Fig. [2fa)), are related by ui = e^^^,r^\ ^.g. 
quired by the Bloch theorem, where K^^n = + liwijX 
and kx G (— tt/A, tt/A). The displacement is fixed to zero 
on the base boundary and the height L of the cell is set 
to 100 /im, two orders of magnitude greater than the sys- 
tem's periodicity. This condition is required to achieve 
a density of states fine enough to appreciate the surface- 
confined system's normal modes also for p <C 1 and p ~ 1, 
the eigenfrequencies spacing scaling as 1/L. Attention is 
devoted in choosing the finesse of the quadratic ordered 
mesh, enabling accurate displacements calculations de- 
spite: (a) the size discrepancy between the nanostruc- 
ture and the silicon cell; (b) the two orders of magnitude 
difference between L and A. Unless otherwise stated, 
the calculations are performed with n = 1 and /c^^ = 
{Kx,\ = 27r/A), corresponding to the first harmonic at 
the center of the surface Brillouin zone, calculation of the 
eigenfrequencies as a function of the wave vector parallel 
to the surface of the substrate not being the scope of the 
paper. 



TABLE I: Material properties for Si substrate^^ and Ni 
stripes:^^ Young's modulus Poisson's ratio a and mass 
density p. 



E (GPa) a p (Kg/m^) 

Si substrate 131 0.27 2330 

Ni stripes 219 0.31 8900 



dj [Cijmn (r) dnUm] = -p (r) u'^Ui 



(2) 
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III. PSEUDO-SAW 

On the contrary to bulk modes, SAWs are confined 
in a tiiin layer starting at the surface of the substrate. 
As seen in Fig. [2fb), in the case of a pure silicon slab, 
without stripes on top, the displacement field vanishes 
within a depth of the order of 1 /im, i.e. of the SAW 
wavelength, corresponding to the chosen grating period 
A. The total displacement color scale normalized over 
its maximum value is reported, together with the dis- 
placement field. These modes are the elastic analogues 
of the electronic surface states^^ or the electromagnetic 
evanescent waves encountered at the surface of a semi- 
infinite solid^ The surface wave is an exact solution of 
the eigenvalue problem, two-fold degenerate under trans- 
lational symmetry. The two solutions u^^ and u^, have, 
respectively, sin {ungerade) and cos (gerade) symmetry. 
In Fig. Efb) only the sin symmetry solution is reported. 
The calculation of the two- fold degenerate eigenfrequency 
at the Brillouin zone center (n = 1) gives z> = 4.92 GHz, 
in agreement with data for a SAW propagating along 
Si(100).24 

On a uniform fiat surface a surface wave does not radi- 
ate into the bulk and has a penetration depth ^ The 
presence of periodic nanostructures is expected to mod- 
ify the properties of the SAW, perturbing the stress and 
velocity fields associated with SAW propagation. ^^^^ In 
the periodic composite system the SAW modes are not 
solutions of the eigenvalue problem. The eigenmodes of 
the full structure, whose displacement field most closely 
resemble the unperturbed SAWs, are shown in Fig. [2fc) 
and Fig. [2fd). The two- fold degeneracy is lifted and a 
band gap opens: one mode, with sin symmetry, is found 
at i>u = 4.27 GHz, whereas a second mode, with cos 
symmetry, is found at = 4.60 GHz. The calculated 
displacement field is Rayleigh-like in region B - evident 
from the alike displacement fields of sin-like solutions in 
region B of Fig. E^b), Fig. Efd) and Fig. Eff) - and bulk- 
like in region C - clear from the field distribution in region 
C of Fig. [2fc) and Fig. [2fd). We address these modes as 
pseudo-SAWs. Periodic nanostructures force the previ- 
ously unperturbed surface waves to radiate elastic energy 
into the bulk. The analogue in the electromagnetic case, 
is diffraction of surface plasmons into far field by a peri- 
odic metallic grating deposited on the substrate surface. 
In the acoustic case, the stress at the Ni-Si interface, 
needed to force the stripes to follow the motion of the 
surface, acts as the source of energy radiation into the 
bulk, as explained by Lin et al.^^ In terms of scatter- 
ing of unperturbed SAWs, this effect can be rationalized 
as coupling of the free surface modes to the silicon bulk 
modes. 

The eigenmodes for the case of full substrate cover- 
age {p = 1) are reported in Fig. [2fe) and Fig. [2ff). 
The band gap closes, the degeneracy is recovered with 
z> = 4.12 GHz, and the eigenvectors are SAWs, hence no 
energy is radiated in the bulk. The solution is similar to 
that reported in Fig. [2fb) for p = 0; the overlay down- 




FIG. 2: (Color online) Nickel stripe on silicon substrate, 
(a) The top 2 /xm portion of the 2D rectangular unit cell 
is reported. The cell geometry is divided in three parts: A. 
Ni stripe, B. 1 /im top portion of Si substrate, C. 99 /im thick 
Si bulk portion. Bloch boundary conditions are set on sides 
1 and 2 of the substrate. The periodicity is A = 1 /xm. The 
height of the Si cell is 100 fim. The Ni stripe is /i = 50 nm high 
and, in this figure, d = 320 nm wide. In Figures (b) trough (f), 
the deformation and arrows correspond to the displacement 
field, while the color scale refers to the normalized total dis- 
placement, (b) Sin-like SAW solution for a pure silicon slab. 
The calculated two-fold degenerate SAW eigenfrequency is 
z> = 4.92 GHz. The orthogonal cos-like solution is not re- 
ported, (c) Cos-like and (d) sin-like pseudo-SAW solutions 
for d = 320 nm. The two pseudo-SAW eigenfrequencies are 
i)g = 4.60 GHz and z>^ = 4.27 GHz. (e) Cos-like and (f) sin- 
like SAW solutions for substrate full coverage. The two-fold 
degenerate calculated SAW eigenfrequency is i> = 4.12 GHz. 



shifts the value of u. 

In the previous discussion, the concept of pseudo-SAW 
was suggested on the basis of similarities between its dis- 
placement field and that of a pure SAW. Nevertheless 
the theory lacks a formal definition for a pseudo-SAW, 
or a procedure allowing to spot out such solutions from 
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TABLE II: Calculated pseudo-SAW frequencies, correspond- 
ing SAW-likeness coefficients and linewidth 7 (FWHM) for 
both gerade and ungerade distributions. 



(a) 



0.4- 
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0.1 


0.32 


0.9 


V, (GHz) 


4.72 


4.60 


4.00 




0.38 


0.12 


0.70 


l9 (GHz) 


0.064 


0.231 


0.013 


vu (GHz) 


4.59 


4.27 


4.15 


Oiu (pu) 


0.18 


0.07 


0.73 


In (GHz) 


0.101 


0.300 


0.010 



tfie infinite eigenvalues set satisfying Eq. [2] in tfie present 
geometry. With this aim in mind, we define the SAW- 
likeness coefficient as 



OL{y) 



{Ea{v)) ^ {EB{y)) 



(3) 



where {Ea{j^))^ {Eb{i^)) and {Etotij^)) are the time- 
averaged mechanical energy contents of u (r) in region A, 
B and in the entire unit cell, respectively (see Fig.[2fa)). 
The calculation is here performed starting from u(r): 



p (r) (v)d^v + p (r) (v)d^v 



(4) 



The SAW-likeness coefficient a (u) outlines which 
eigenmodes of the system have mechanical energy 
mainly localized within a depth A (where A is both 
the nanostructure's period and the penetration depth 
of an unperturbed SAW of the same wavelength). In 
Fig. [31 we report a {u) for three values of p. Two 
distributions, ag {u) and (i^), emerge from the plot, 
corresponding to the cos-like and sin-like displacement 
profiles, respectively. We define pseudo-SAWs the 
eigenmodes with frequencies i>g and z>n, corresponding 
to the two distributions' maxima. The pseudo-SAWs are 
the solutions that most likely resemble the unperturbed 
SAWs. The eigenfrequency degeneracy is removed 
and the opening of a frequency band gap is directly 
shown. We pin-point that the SAW-likeness coefficient is 
easily interpreted as the SAW line-shape, the line-shape 
concept arising because, in the composite system, the 
SAW is not an eigenmode. 

In Fig. [3fa) we report the calculation for a sample 
with p = 0.1 {d = 100 nm). For small filling fractions, 
nickel stripes act as a weak perturbation to the modes 
of the underlying substrate. The values of the pseudo- 
SAW frequencies and the corresponding SAW-likeness 
coefficients, as well as the linewidth 7 (FWHM) of the 
two distributions, are reported in Table [III The fact 
that OLg{vg) > au{i^u) means that the grating affects the 
sin-like SAW more than the cos-like one. Furthermore, 
the grating couples the unperturbed sin-like SAW over 
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FIG. 3: (Color online) SAW-likeness coefficient a (jy) ver- 
sus the calculated eigenfrequencies. Two distributions arise, 
au (i^) and ag (zy), corresponding to sin- like (empty circles) 
and cos-like (crossed circles) displacement profiles. Calcula- 
tion for configurations with different filling fraction are re- 
ported: (a) p = 0.1, (b) p = 0.32 and (c) p = 0.9. Solid lines 
are intended as a guide to the eye. 



a wider range of modes: 7^ > 7^. The frequency 
gap is Az> = i>g — = 0.13 GHz. Both Dg and Du 
are lower than the unperturbed SAW eigenfrequency 
z> = 4.92 GHz. 

Increasing the filling fraction to p = 0.32, we 
obtain the distributions reported in Fig. ^h). The 
physics is now well beyond the perturbative approach: 
larger linewidths and lower values of the SAW-likeness 
coefficient for the two pseudo-SAWs are the signatures 
of a stronger SAW-stripes interaction, as compared to 
the case of p = 0.1. This leads to a higher coupling 
of surface waves with the bulk modes of the system, 
a confirmation to the evidence given in Fig. [2fc) and 
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Fig.[2];d). The frequency gap is Ai) = 0.33 GHz. Both Ug 
and i>u are further down-shifted with respect to z>. The 
pseudo-SAW frequency lowering is strongly dependent 
on the displacement field symmetry. 

Increasing the filling fraction to p = 0.9, close to 
silicon full coverage (see Fig. [3fc)), narrow linewidths 
and high peaked distributions resemble the frequency 
two-fold degenerate surface wave solutions shown for 
p = 1 in Fig. ^e) and Fig. [2ff), consistently with a 
reduced frequency gap Az> = —0.15 GHz. Interestingly 
Az> changes sign and the ungerade pseudo-SAW is now 
less affected by the periodicity as compared to its gerade 
counterpart. For the case of p = 1, z> is in agreement 
with the analytical theory developed by Auld.^^ 

To further test the soundness of the pseudo-SAW 
definition adopted in the present work, we calculated for 
p= 0, 0.32 and 1, the dispersion for the first harmonic 
pseudo-SAW modes: e [-0.27r/A, 0.27r/A] and 

n = 1. The outcome is shown in Fig. |4j The dispersion 
relation for the semi-infinite silicon slab {p = 0) starts at 
u = 4.92 GHz and is linear with a slope Vg = 4900 m/s, 
in agreement with the speed of sound reported for SAW 
on Si(lOO).^^ For the fuh overlay case (p = 1), the dis- 
persion starts at 4.12 GHz. Both solutions are true SAW 
eigenmodes of the system. For the case of p = 0.32, we 
observe a gap opening and the two dispersion branches 
for the pseudo-SAWs remain correctly confined between 
their p = and p = 1 analogues. 



IV. FREQUENCY GAP 

Having formally defined the pseudo-SAW, we now ad- 
dress the issue of the frequency gap in the pseudo-SAW 
modes. In Fig. [5] the pseudo-SAW frequencies and band 
gap calculated over the entire filling fraction range are 
shown. For p 7^ 0, the pure silicon slab eigenfrequency 
degeneracy is lifted. As p increases, Dg decreases un- 
til it reaches a plateau for p = 0.2. The plateau lasts 
for p G (0.2, 0.5) then, for p > 0.5, Ug decreases mono- 
tonically till it jumps back to the SAW value for full 
coverage as soon as the stripes' width d reaches the size 
of the periodicity, becoming a complete overlay. There 
is a geometrical explanation for the discontinuity (see 
inset of Fig. [5|). For p = 1, periodic boundary con- 
ditions become effective also on the stripe's sides, the 
cos-like pseudo-SAW solution undergoes a transition 
from free to constrained displacement on the stripe's lat- 
eral boundaries. The same constraints do not affect u^^, 
because they are set where the nodes of the sin-like solu- 
tion are found (see Fig. [2l^f)). In an experiment, for the 
above mentioned discontinuity, a steep positive deriva- 
tive dvg/dp for p ^ 1 should be expected in place of an 
abrupt transition. Concerning v^-, it is a monotonically 
decreasing function of p. It decreases with a steep slope 
for p < 0.3. It then changes slope stabilizing to an al- 
most constant value all the way io p = 1. For p < 0.75, 
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FIG. 4: (Color online) Dispersion relations of the elastic sur- 
face modes for small values of kx (in tt/A units) and n = 1, 
for a pure silicon slab (p = 0), for silicon full coverage (p = 1) 
and for the configuration with Ni stripes over Si substrate 
(p = 0.32). Solid lines are fits to data. Dotted lines are a 
guide to the eye. 



> at p = 0.75 the eigenvalues cross; for p > 0.75, 
i>g < Vu] degeneracy is recovered for p = 1. The fre- 
quency gap Az> opens as soon as p ^ 0. The maximum 
value Az> = 0.40 GHz is attained when the silicon surface 
is close to half coverage. For higher filling fractions, the 
gap decreases to Az> = for p = 0.75. For values of p 
in excess of 0.75, Az> < and decreases till it makes an 
abrupt transition back to Az> = for p = 1. The present 
results also show the limits of applicability of the per- 
turbative approach. Within a perturbative scheme the 
relative frequency shift is proportional to the filling fac- 
tor: (z^(p) — i>(0))/z>(0) oc {h/X)p. Calculation of the rel- 
ative frequency shift, on the basis of the results reported 
in Fig. O shows that linearity holds up to p ~ 0.1. For 
values of p in excess of 0.1 one has to rely on the full 
calculations here reported. 

In the general case, for a fixed periodicity A, 
the pseudo-SAW frequencies - both gerade and unger- 
ade - are functions of the parameters {p, /i,m,p}: z> = 
f {p^h^m^ p). The four parameters are not mutually in- 
dependent, being p oc m/{hp). v is then defined in a 
four-dimensional space and it is a function of three pa- 
rameters only, the fourth being fixed by the choice of the 
other three. For instance, each data point reported in 
Fig. [5] (left axis) is calculated having fixed p, h and p 
independently. The plots are then obtained spanning p 
over the range (0, 1), thus exploring a particular trajec- 
tory z>(p, = 50 nm, pNi) on the hyper-surface z>(p, /i, p). 
A full description is computationally too time consum- 
ing and would not add much to the comprehension of the 
physics without slicing the hyper-surface along particular 
hyper-planes or trajectories. 

To get a deeper physical inside, the effects on z> of 
the geometric factors {p, h} and mass loading {m, p}, are 
here explored. We start considering a reference system 
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filling fraction p (or stripes width d in nm, top axis). Right 
axis: pseudo-SAW frequency gap Az> (triangles). The hori- 
zontal line is the Az> = line and refers to right axis only. 
Inset: magnification of the graph for high filling fractions p is 
shown. Lines are a guide to the eye only. 



I 



0.6-' 



0.4- 



0.2- 



(b) 



o 



\_ 



0- 



0.5- 
0- 
-0.5- 
-1.0- 
-1.5- 
-2.0- 







c 

L 


) COS symm 
) sin symm 
^ frequency 


etry 
etry 
gap 


> ^ 


5 e 


5 


-4 





\ J< 












s) 










5 <k 


^ ^ 


5 ^ 


0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Filling fraction p 



0.2 



0.1 CD 

c 

CD 
CQ 

-0.1^ 

o 

-0.2 



--0.3 



® COS symmetry 
o sin symmetry 
A frequency gap 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Relative stripes mass density p/p^^^ 



configuration with p = 0.2, = 50 nm and p = pNi-, 
hence the value of m is fixed to m^^^^\ The pseudo- 
SAW frequencies, calculated in the reference configura- 
tion, are i/^^^*^^ and Vu^^^ • To test the contribution of 
the stripes' geometry, m and p are kept fixed at their 
reference values, while p is allowed to increase from 0.2 
to 0.9. Doing so, the value h decreases. In Fig. [6l^a) 
the dependences of i>i — i)^^^-^^ - with i = {g^u} - and 
Az> on the geometric factor {p, h} are shown. The geo- 
metric factor up-shifts both frequencies, Uu being more 
affected than i>g and Az> decreasing monotonically. The 
effect of the geometry in surface-based phononic crystal 
fits well with the results reported by Tanaka et al^^ for a 
2D phononic crystal in which no mass loading effect is to 
be expected (the band gap definition is reversed in sign 
with respect to our case). This makes us confident that 
we are properly disentangling the geometry from mass 
loading effects. To inspect the mass loading effect, p and 
h are kept fixed at their reference values, while p is raised 
up to 4.5 times the Ni density. Doing so, the value of m 
increases. In Fig. [6fb) the dependences of i>i — i)^^^^^ 
and Az> on the mass loading factor {m, p} are shown. 
The mass loading down-shifts both frequencies, Vu being 
more affected than Vg and Az> increasing monotonically. 
The functions in Fig. [6l^a) cannot be compared with the 
ones reported in Fig. [6fb), nor the absolute values are 
of relevance. The importance stands in the qualitative 
functions behavior, showing the disentangled physical ef- 
fects of geometry and mass loading over z>^ and Az>. The 
present analysis shows that in a device the pseudo-SAW 
frequency gap results from a combination of geometry 



FIG. 6: (Color online) Geometric and mass loading effects 
on z>. The reference system configuration is p = 0.2, /i = 50 
nm and p — pNi, hence m is fixed to inrS^^^\ The pseudo- 
SAW frequencies, calculated in the reference configuration, 
are Vg^^^^ and i)u^^^ ■ (a) Exploring the geometric effect: the 
mass loading factor {m, p} is kept at the reference value, the 
variable being the geometric factor {p, /i}. (b) Exploring the 
mass loading effect: the geometric factor {p, h) is kept at the 
reference value, the variable being the mass loading factor 
{m,p}. Left axis: pseudo-SAW eigenfrequencies Vg (crossed 
circles) and i)u (empty circles). Right axis: pseudo-SAW fre- 
quency gap Az> (triangles). The horizontal line is the Az> = 
line and refers to right axis only. Lines are a guide to the eye 
only. 



and mass loading factors, the two affecting the pseudo- 
SAW frequencies in opposite ways. A trade-off between 
p and m has to be foreseen for tailoring the pseudo-SAW 
frequencies or frequency gap in view of possible applica- 
tions. 



V. MECHANICAL ENERGY RADIATION 

In this section, the energy distribution in the system, 
for both pseudo-SAW eigenmodes and u^, is studied 
as a function of p. Energy distribution in surface-based 
phononic crystal can be attained from the knowledge of 
the displacement or velocity and stress fields. The fields 
are typically calculated via coupled-mode theory or scat- 
tering theory in Born approximation. These are pertur- 
bative approaches. For instance, in this latter approach 
the SAW solution for the half-infinite Si slab is taken 
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as the zero-order wave (impinging wave) used to eval- 
uate the scattering matrix - the interaction term being 
brought in by the grating - together with final states 
modes (scattered waves). Such an approach loses relia- 
bility as the SAW frequency increases. Since the pene- 
tration depth of the SAW equals its wavelength A, higher 
SAW frequencies imply stronger surface confinement and, 
ultimately, stronger scattering with the periodic grat- 
ing. This requires considering higher order terms in the 
scattering integral equation in order to achieve a reason- 
ably precise solution. For surface-based phononic crys- 
tals working in the hypersonic range (A < 1 /im) the 
strong confinement is an issue. The theoretical frame 
here introduced to define the pseudo-SAW, together with 
the calculations of the composite system eigenmodes, can 
be exploited to tackle the problem of energy distribution 
in the system beyond a perturbative approach and the 
results can be readily translated in a scattering frame- 
work. 

The definition given in Eq. [3] is here extended to ac- 
count for the time-averaged normalized energy content 
in the three regions outlined in Fig. [2]^a) and in the inset 

of Fig. m 

^-(^) = 7^4%' x = {AB,C} . (5) 

In Fig. [3 the energy contents ax{i>g) and ax{i>u) are 
reported for x = {A, 5,(7} as a function of the filling 
fraction p. The results show the mechanical energy spa- 
tial distribution - normalized against the total energy in- 
side the cell - over the entire filling fraction range, when 
a pseudo-SAW eigenmode is excited. For p = the SAW 
is a true eigenmode, the energy being concentrated, as 
expected for a SAW, in region B: a^(z>) ~ 1. For small 
filling fractions {p < 0.05), at least for u^, the elastic 
energy is mostly concentrated in the 1 /im top portion 
of the Si slab, where the unperturbed SAW is expected 
to dump. The grating is much more effective in scat- 
tering energy to the bulk for u^^: ^^(^n) ^ <^c(^n)- 
This symmetry-related difference in coupling a SAW to 
the bulk is a confirmation to the observed relationship 
^ of Fig. [3l^a). As p increases to half cover- 
age, a greater amount of mechanical energy is transferred 
from the Si surface region B to the bulk region C. From 
the eigenvalue problem perspective, a greater amount of 
energy radiated into the bulk translates in energy dis- 
tributed over a wider range of eigenmodes of decreas- 
ing SAW- likeness coefficient a(z>), 7 increasing with p for 
both gerade and ungerade pseudo-SAWs (see Fig. [3l^b)). 
The maximum values for both ac{i^g) and ac{i^u) are 
attained for p = 0.4, the same filling fraction maximizing 
Az>. For p > 0.4 the trend reverses, being dac{i^)/dp < 
whereas daA,B{i^)/dp > 0; the relative energy content in 
the bulk region is transferred in region B and A. For 
p ^ 0.8 the energy content in the Ni stripe equals the 
energy content in region B; the equality of energy con- 
tent holds also between sin and cos symmetry eigenmodes 
and pseudo-SAW frequency crossing occurs (see Fig. [5|). 
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200 400 600 800 1000 
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FIG. 7: (Color online) Analysis of the spatial distribution 
of the mechanical energy in the composite system over dif- 
ferent filling fraction configurations, normalized against the 
total energy inside the cell. In the inset, the circles represent 
the fraction of the mechanical energy confined in the nickel 
stripes region A, the triangles are for the top 1 pm portion B 
of Si substrate and the squares are associated with the silicon 
bulk C. The dotted lines are a guide to the eye only. 

As the Si full coverage configuration is approached, the 
radiation into the bulk is further reduced in favor of a 
strong mechanical energy confinement in the stripe and 
in the top 1 jam portion of Si substrate, the crossing point 
being for p = 0.9. For p > 0.9, the oscillation of the 
Ni stripe starts to energetically dominate and the stripe 
behaves as a low- loss acoustic waveguide. The unper- 
turbed SAWs, eigenmodes of the semi-infinite slab, are 
now heavily scattered into the Ni overlay. In the limit of 
Si full coverage (p = 1), the Ni overlay acts as a lossless 
waveguide. The SAW is a true solution of the eigenvalue 
problem, without coupling with bulk modes, and it is 
frequency down-shifted with respect to the overlay-free 
SAW by an amount ~ 0.8 GHz. 

VI. CONCLUSIONS AND PERSPECTIVES 

In the present work we propose a theoretical frame 
allowing to access the physics of pseudo-surface acoustic 
waves in surface phononic crystals. The pursued strategy 
can be applied to any surface phononic crystal, enabling 
investigation of pseudo-SAW line-shapes, gap opening 
and mechanical energy scattering beyond perturbative 
approaches, thus finding application also in the hyper- 
sonic frequency range. 

We applied the outlined theoretical framework to the 
case of a hypersonic surface phononic crystals made of 
periodic nickel stripe on a silicon substrate. We then 
investigated the pseudo-surface acoustic wave frequency 
gap over the entire filling fraction range, starting from 
the case where the periodic nanostructures act as a per- 
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turbation, to finally reach the substrate full coverage. An 
understanding of how the construction parameters affect 
the frequency gap has been achieved. We showed that the 
pseudo-SAW frequency gap results from a combination of 
geometry {p, h} and mass loading {m, p} factors, the two 
affecting the pseudo-SAW frequencies in opposite ways. 
A trade-off between p and m is thus necessary for tai- 
loring the pseudo-SAW frequencies or frequency gap for 
the application at hand. The mechanical energy spatial 
distribution of pseudo-SAW as a function of the filling 
fraction has been studied, allowing to tailor the device's 
parameters in order to minimize the energy content scat- 
tered out of the desired modes. 

The outlined theoretical scheme will proof a useful 
tool for future applications involving hypersonic pseudo- 
SAW in surface-phononic crystals. Complete control 
of the phonon gap opening in hypersonic surface-based 
phononic crystals, addressed in the present work, will be 



of impact in various fields of modern solid state physics. 
For instance the acoustic gap can be exploited to engi- 
neer surface waveguides and cavities to guide and localize 
hypersonic pseudo-SAWs. In addition, application of the 
recently demonstrated negative refraction effect^^ to sur- 
face phononic crystals, will foster the design of acoustic 
lenses working at hypersonic frequency. The possibil- 
ity to focus and localize hypersonic pseudo-SAWs will 
enable efficient manipulation of the optical, electronic 
and magnetic properties of nanostructured devices in lo- 
calized spatial regions. Finally, recent advances in the 
near-field optical microscopy pave the way to the cou- 
pling of the time-resolution of standard pump-probe op- 
tical experiments^^^^^ to the spatial resolution of SNOM 
probes. This will allow to study in real-time the local- 
ization of hypersonic pseudo-SAWs in nanostructures or 
surface cavities with sub- wavelength spatial resolution. 
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